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We give a summary of the theoretical attempts to make sense of the observational data concerning the highest 
energy cosmic rays, i? > 3 x 10^^ eV. 



1. Preamble 

It was about forty years ago that Linsley PP 
presented an incredulous physics community with 
the ^ 10^" eV event observed at the Volcano 
Ranch array. Today ground arrays and fluores- 
cence detectors have collected over a hundred 
events with energies above 3 x 10^^ eV, and a 
dozen or so at or above 10^*^ eV . It is worth 
remembering that the ~ 2 TeV CM energies 
attained at the Tevatron are at the knee of the 
cosmic ray spectrum. LHC will operate at an 
equivalent fixed target energy falling approxi- 
mately halfway between the knee and the ankle, 
around the so-called second kne^^ region. Thus, 
above the second knee, the first interaction in an 
extensive air shower involves physics untested by 
the well-controlled experimental environments of 
accelerators. 

For many years, the appearance of the GZK 
feature [21 at the highest energies (i.e. above 
2 X 10^^ eV) was taken for granted. The pro- 
cesses, p + 7cMB + X, and, of lesser impor- 
tance, p + 7CMB e+e^ -f X, are well-studied 
low energy processes in the CM. A proton of 
energy 10^° eV or above has a mean free path 
of 7-8 Mpc in the cosmic microwave background 
(CMB) radiation. With each interaction event, 
the proton loses 20-25% of its energy; after only 
a few interactions its energy drops under the 
threshold of the pion photoproduction process. 
The expected precipitous drop in the incoming 
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cosmic ray flux at 5 x lO^^ - 8 x 10^^ eV (GZK 
cutoff) and the pileup |3] (a bump in the spec- 
trum) just before the drop"^ are called "the GZK 
feature" . The GZK sphere approximately defines 
the cosmic neighborhood of sources from which 
we expect the arrival of extreme high energy 
cosmic rays [E w 10^° eV and up). We usually 
put this particle horizon at 50-70 Mpc. However, 
the GZK radius could be considerably less as 
demonstrated by Stanev et al. H. Scattering 
in extragalactic random magnetic fields makes 
the real distance travelled by the particle longer 
than the source-Earth distance; more energy loss 
processes can take place as a result. Heavier 
nuclei with similar energies also suffer energy 
loss through photodisintegration on the cosmic 
microwave and on the intergalactic infrared back- 
ground radiations For ^^Fe, for example, the 
decrease of the energy loss time is even more 
abrupt than for protons, though it starts at a 
somewhat higher energy. 

We will summarize the most important, and 
often controversial, observational data. Then we 
present the theories formulated to explain certain 
coherent subsets of the data. 

2. Observational Data 

What follows is a very brief summary of the 
most important features of the extensive air 
shower and EGRET data from a theorist's point 
of view. Only ultra high energy cosmic rays 
(UHECR) with energies above 3 x 10^^ eV are 

■^This is a consequence of the pion photoproduction by the 
above threshold protons or nuclei. 
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considered. 

2.1. Spectrum 

At and above 10^° eV, there are 15 events: 11 
observed by AGASA and 2 by HiRes |7| and 
one each by Fly's Eye and Yakutsk 

There is a substantial disagreement between 
AGASA and HiRes at this time. The HiRes data 
were successfully fitted with a two-component 
model jl(J| including primaries with galactic and 
extragalactic origins. The Fly's Eye data sug- 
gest a transition from heavy to light (proton) 
primaries (at energies from 10^^ to 10^^ eV), 
which can be interpreted as a sign of the passage 
from galactic to extragalactic sources. The extra- 
galactic source model of Berezinsky et al. |12| . 
which takes into account both uniform source 
distribution and a maximal energy at 10^^ eV, 
yields not only the GZK feature but also other 
details of the spectrum, as a result of considering 
both energy loss processes. The large number of 
events above the GZK cutoff (mostly observed 
by AGASA) remains unexplained by the two- 
component models'*. 

It should be emphasized that even the two 
largest experiments, AGASA and HiRes, have too 
few events at the highest energies; the discrep- 
ancy between them has low statistical significance 
as a result [T^ . 

2.2. Composition 

As mentioned above, fluorescence detectors 
provide some hints for a Fe p transition. How- 
ever, it is important to note that the mass compo- 
sition results are different for each of the various 
methods and model calculations used^l]. In ad- 
dition, the shower energy estimates are sensitive 
to the assumption about the nature of the pri- 
mary. To quote A. A. Watson {loc. cit.), "knowl- 
edge about the mass of the primary cosmic rays 
at energies above 10*^ eV is rudimentary." As 
far as photon primaries are concerned, it is very 

* Since for both AGASA and HiRcs there are systematic 
errors in the energy assignment, it is difficult to judge how 
great the disagreement between them actually is. 



unlikely that the primary of the "gold-plated" 

3 X 10^" eV Fly's Eye event was a photon [T3] . 
Around 10^® eV, the photon flux is less than 
40% of the baryonic component according to both 
Haverah Park JJ, and AGASA .17^ . By search- 
ing for deeply penetrating horizontal air show- 
ers, AGASA investigated the UHE neutrino flux 
for energies > 3 x 10*^ eV ^H]- No significant 
enhancement above the expected hadron-induced 
flux was reported. 

2.3. Directionality 

For the energy range we are interested in, 
the arrival direction distribution is remarkably 
isotropic. This could indicate a large number 
of weak or distant sources. AGASA data show 
small-scale clustering (4 doublets and one 
triplet within a separation angle of 2.5°) above 

4 X lO^^eV (57 events). HiRes in stereo mode 
above lO^^eV (162 events) claims to see no small- 
scale anistropy in the arrival direction distribu- 
tion 1201 • Recently, it was found [23 that the 
statistical significance here is also very low for 
both the small scale clustering in the AGASA 
data set and the apparent discrepancy between 
AGASA and HiRes. 

Nevertheless, despite the paucity of data 
points, statistical analyses have attempted to 
discover possible angular correlations between 
the arrival directions of the extremely high en- 
ergy cosmic rays (EHECR), and to match these 
in turn with high energy source canditates. The 
enticement to carry out an analysis of even the 
small available sample is obvious. Significant 
angular correlation in the arrival directions of 
the EHECR could be the a smoking gun for the 
existence of discrete compact (repeating) sources. 
Using the AGASA data of S > 4 x lO^^ eV ^ the 
presence of clustering can be used to constrain 
the number of sources in the nearby Universe 
for proton primaries p2] to be 10~^Mpc~^ with 
large uncertainties. Tinyakov and Tkachev 
claimed to have identified BL Lacertae objects 
as sources generating the anisotropies in the 

^At E = 4: X lO^^eV protons have a loss length of ~ 
lOOOMpc. 
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AGASA data. Evans et al. ^ and W. Burgett 
et al. questioned this conclusion, however^. 

Most authors working in the field agree that 
more data from HiRes and from larger experi- 
ments like the Southern (and hopefully the North- 
ern) Pierre Auger Observatory, EUSO, the Tele- 
scope Array, and OWL are necessary to draw 
solid conclusions about clustering and the pres- 
ence or absence of the GZK cutoff discrepancies. 
Until then the BL Lac alignment is merely sug- 
gested by the need for objects accelerating par- 
ticles to exceptional energies, rather than being 
a firmly established consequence of the observa- 
tional data. 

2.4. Extragalactic Diffuse Gamma Rays 

Another piece of important data is the ex- 
tragalactic diffuse 7 ray background (EGRB). 
It is generally believed that unresolved AGN, 
galaxy clusters, distant gamma ray bursts, etc. 
are the main sources. However, EGRB would 
get a significant contribution from the decay of 
topological objects, relic particles etc., which are 
used by some authors to explain trans-GZK cos- 
mic ray events. Comparing that contribution to 
the measured EGRB, one may be able to place 
upper limits on, or eventually exclude, some of 
those schemes. Any electromagnetic (EM) en- 
ergy introduced to the EM background radiation 
(dominated by CMB) produces an EM cascade 
if the energy is > 10^^eV/(l -I- z) i.e. the pair 
production "threshold" on the CMB at redshift 
z. The cascade by pair production and inverse 
Compton scattering degenerates the injected en- 
ergy to under ~ 10^^ eV in a couple of Mpc and 
produces a power-law distribution of GeV-range 
gamma rays. 

To obtain the EGRB, one must construct a 
model of the galactic 7 ray background, subtract 



^Burgett and O'Malley partition the AGASA sample into 
three statistically uncorrelated groups by energy: E < 
5 X 10^^ eV, 5 X 10^^ < E < S x 10^^ eV, and E > 
8x 10^^ eV. They find the middle group to be preferentially 
aligned with the Supergalactic equatorial plane, while the 
other two groups are statistically consistent with isotropic 
distributions. 



it from the observational data, and assign the re- 
mainder to the EGRB. This is a complicated and 
model-dependent process. The EGRB was first 
discovered with SAS-II and then confirmed 
by EGRET ^7\. A power law energy spectrum 
(^(-(2.i±o.o3)) ^-j^g j.j^^gg 30 MeV < E <100 

GeV) with no anisotropy was found. Strong et 
al. |2H|, in a recent reworking of the EGRET 
data, substantially improved the determination 
of the diffuse galactic 7 ray flux. They conclude 
that in previous analyses the galactic contribu- 
tion was underestimated; the EGRB must be 
revised downward as a result. The new spectrum 
is not consistent with a power law. One may infer 
tentatively that the EGRB arises from a number 
of different physical processes. 

At higher energies (20TeV < < SOOTeV), 
an upper limit on the EGRB can be obtained by 
estimating the ratio of photon-induced to hadron- 
induced extensive air showers. The GRAPES 
collaboration |29| obtained an upper limit of 
about 3 X 10"'' for the fraction of photon-induced 
(muon-poor) air showers. 

3. Theoretical Models 

The presence or absence of the GZK feature 
determines whether it is enough to rely on con- 
ventional astrophysics to explain the spectrum, 
composition, and directionality encrypted in the 
data, or if physics beyond the Standard Model is 
needed IHHj . 

3.1. Traditional Astrophysics 

These models are also called bottom-up scenar- 
ios. At 10^^ eV and above, the extragalactic 
component must dominate the energy spectrum. 
Even if there are some objects (neutron stars, 
galactic wind, etc) in our galaxy which could pos- 
sibly accelerate UHECR, it would be impossi- 
ble for them to produce the observed isotropic 
distribution. For order of magnitude estimates, 
the community generally uses Hillas' formula. 
Whether the acceleration takes place electromag- 
netically or with the help of diffusive shock accel- 
eration (first order Fermi process), the magnetic 
field must be strong enough to keep the particles 
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in the region of acceleration, thus, approximately, 

R > R ~ ^ 

-'I'accclcrator ^ -^gyromagnctic ^ rj t-j i 

where the charge of the accelerated particle is Ze. 
From this. 

We are looking for a maximal energy of about 
10^^ — 10^^ eV. The enormous energy of the 
primary singles out exceptional astrophysical 
sources: radio lobes of powerful radio galaxies, 
jets of active galactic nuclei, gamma ray bursters, 
and magnetars. There are only a few of these 
possible accleration sites inside the GZK radius 
for charged hadrons. (For example, the radio 
galaxy M87, in the center of the Virgo cluster, 
was investigated repeatedly as a candidatefST]. 
It could provide EHECR at Earth only if the ex- 
tragalactic magnetic fields on the way here have 
some specific configuration.) Consequently, it is 
unlikely that the GZK cutoff could be completely 
absent in the spectrum. The accelerators are 
far beyond the pion photoproduction distance 
and the random, poorly known, extragalactic 
magnetic fields may cause deflection of charged 
particles. In fact, the relevant sources are at cos- 
mological distances. Nevertheless there are small 
scale anisotropics and correlations with certain 
sources. These observations indicate that the 
cosmic traveler may be a stable, neutral parti- 
cle which interacts very weakly with the cosmic 
background radiation fields. This particle could 
be a neutrino. 

It is generally accepted that the dominant form 
of acceleration of cosmic rays at the sources is 
by first order Fermi acceleration at astrophysical 
shocks. A recent development in the theory of 
shock acceleration is the introduction of the back 
reaction of the accelerated particles to the shock 
modifying it stronglv|32j. Among other results 
the spectrum is claimed to be flatter than in the 
linear theory; thus, the number of arriving trans- 
GZK primaries could be higher. 



3.2. Scenarios Beyond the Standard 
Model 

3.2.1. Z-bursts 

The Z-burst modeljnHl needs only the inclu- 
sion of the neutrino mass in the Standard Model 
(SM) for its attempt to explain the trans-GZK 
events. Among SM particles only the neutrino 
can handle cosmic distances without dramatic ab- 
sorption. The assertion is that extreme energy 
neutrinos 

= 4 X {m^{eV))-^ZeV 

forming bosons resonantly with antineutrinos 
of the relic 1.9 °K thermal neutrino background 
could produce enough protons well inside the 
GZK radius to account for the EEGR. On aver- 
age a produces ~ 2 nucleons, ~ 20 photons 
and ~ 50 neutrinos. Although the Z-burst model 
is very economical, but it is likely to have fatal 
problems both because of the required energy 
at the source and because of the extreme en- 
ergy neutrino flux required to produce enough 
of the trans-GZK events PIJ. The new experi- 
mental hmit by WMAP, 2dF and SDSS on the 
active neutrino mass (^m^i < 0.72 eV) com- 
bined with the neutrino oscillation results gives 
the maximal possible neutrino mass rrii, « 0.24 
eV [33]. This puts the energy of the neutrino for 
the resonant production at ~ 10 ZeV. If these 
neutrinos are secondaries of pions, the protons in 
the source must be accelerated to ~ 10^ ZeV. No 
known astrophysical object is capable of acceler- 
ating protons to such high energies [21] • A large 
enough flux by ordinary astrophysical sources (or 
by beyond-the Standard Model objects [23) to 
produce a substantial part of trans-GZK events 
would pump so much EM energy into the uni- 
verse that the EGRET limits would be violated 
many times over. 

An overdensity of neutrinos in a halo of our 
galaxy or in that of the Local Group could help 
to reduce the flux requirement. However, these 
neutrino masses are too small to allow any clus- 
tering of relevance to allow the reduction of the 
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incoming neutrino flux [SHI- For a conserva- 
tive maximal neutrino mass, rrii, « 0.33 eV, the 
neutrino flux required by the Z-burst scenario 
is aheady in marginal conflict with the upper 
limit given by the FORTE ,^0, experiment and 
in violation of the GLUE |1J limit. A lower 
mass makes the problem more severe |3y| . This 
limit on the Z-burst model is valid even when the 
distant sources emit only neutrinos. 

3.2.2. Strongly Interacting Neutrinos 

The problem of energy loss during the propa- 
gation at cosmological distances would not arise 
if the trans-GZK primaries were neutrinos and 
interacted strongly at very high energies |42| . 
causing proton- like air showers. Both the energy 
and flux reqirements are much lower here than 
in the Z-burst model. In this case the cosmo- 
genic neutrino flux (i.e. secondaries from the 
pion photoproduction of extreme energy protons 
on the CMB) could account for the trans-GZK 
events |43| . since all incoming neutrinos produce 
an air shower. For the protons the injection spec- 
trum is simple, power law. The onset of strong 
interaction must result in a very rapidly rising 
cross section; a slow rise would give a large num- 
ber of horizontal showers, which are not seen by 
either AGASA or HiRes. For example, in the- 
ories of higher than 4 dimensions, Kaluza-Klein 
graviton exchange gives a power law increase, 
which is not only too slow, but strong interaction 
cross sections cannot be reached until energies 
« 1020 eV EH. 

The following two schemes also assume some 
higher dimensional (string inspired) theory. The 
neutrino-nuclcon cross section begins to grow at 
a certain energy, in essence due to the unification 
of all interactions at a scale of the order of a few 
tens of TeV or a few TeV in the CMS, producing 
lepto-quark resonances 03] , mini black holes 0^1 
or "p-branes" 07j. The asymptotic (s ^ oo) 
equivalence of the resonance and mini black hole 
schemes was recently demonstrated EH]- The 
onset of the new interaction is very different. 
In the "semiclassical black hole" picture, it is 
assumed that the cross section is proportional 



to the square of the higher dimensional (d=10 
?) Schwarzschild radius. Combined with the 
assumption that mbh oc ^/s, it leads to an un- 
acceptably slow rise of the cross section. The 
published data on near horizontal showers by 
AGASA and Fly's Eye exclude that. By con- 
trast, with a string scale of the order 70 TeV one 
is far below the semiclassical region jj^]- The 
growth of the cross section is governed by the 
exponentially increasing level density in the tran- 
sition region (and only in the transition region) 
from the SM to the string regime. In this scheme 
few near horizontal showers are produced, as a 
result. 

Another possibility, suggested to achieve strong 
cross sections, is the inclusion of instanton in- 
duced processes. There is a detailed discussion 
and further references on the instanton enhanced 
cross section in |43| . 

It is remarkable that the cosmogenic neutrino 
beam, together with a strongly interacting neu- 
trino with sudden onset of large cross section 
could give a solution of the trans-GZK events. 
This is accomplished without the need to as- 
sume exotic astrophysical scenarios. There is no 
obvious contradiction with data obtained from 
FORTE and GLUE, either. 

3.2.3. Top-down models 

Very heavy topological defects (TD) could 
solve both the problem of the long trip for pro- 
tons in the CMB and the difficulty to accelerate 
particles to extreme energies [HI]. In the early 
universe at the time of phase transition from the 
GUT scale (M > lO" GeV) TD (monopoles, 
cosmic strings, necklaces, textures etc.) could be 
formed in localized regions where mass-energy is 
trapped. In the decay (or annihilation) of these 
superheavy extragalactic relics, "X-particles" are 
formed which decay trough QCD fragmentation 
jets consisting mostly of photons, neutrinos and 
~ 5% nucleons. These nucleons provide the 
trans-GZK events. However, a large fraction of 
the energy is injected into the extragalactic space 
via 7-s with E > 10^^ eV. Thus, if the TD can 
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generate the flux of EECR, then the EGRET 
hmit (see 2.4) is always exceeded jSU) . 

Another type of top-down scheme uses the 
possible existence of superheavy dark matter 
particles, which could have been created at the 
end of inflation in the early Universe. Their decay 
or annihilition could also produce the trans-GZK 
spectrum via QCD fragmentation As cold 

dark matter, they could have an overdensity in 
the halo of our galaxy by a factor of 10'*. Their 
mass must be > 10*^ GeV and the extreme energy 
protons will reach the detectors since they are 
produced closeby. However, there is anisotropy 
in the arrival direction because of our position in 
the galaxy. In addition, the fragmentation pro- 
cess produces about a factor of 2 more ultra high 
energy photons than protons. As it was discussed 
in Section 2.2, the experimental photon fraction 
is certainly < 1. 

3.2.4. Stable Light Neutral Particles 

The idea of light, supersymmetric, stable or 
long-lived hadrons instead of nucleons as UHECR 
was entertained for some time. A mass larger 
than the mass of the nuclcon (« 2 GeV) would 
increase the threshold energy of the tt photo- 
production and push the GZK cutoff to higher 
energies. 

One hypothetical particle in this class re- 
quired the existence of a light gluino. Precision 
measurements at LEP and SLC closed the so- 
called "light gluino window" ^5.3,^ . Shadrons with 
a light sbottom were investigated as well |54j . 
motivated by the preliminary result from the 
Tevatron indicating an overproduction of bb. 
Futher data taking, however, did not confirm the 
effect ISg. 

The possibility of axion-photon mixing was 
considered in Ref. The hypothesis was that 
the trans-GZK showers were started by regener- 
ated photons, while the axions were the long dis- 
tance travelers. One may object, however, that 
according to present observational data the trans- 
GZK showers look hadronic and not photon- 



induced. 

3.2.5. Lorentz Invariance Violation 

The idea of small violation of Lorentz invari- 
ance is not new (57) . The idea was reconsidered 
again as a possible explanation for the lack of 
the GZK cutoff in the AGASA and Fly's Eye 
data (58) as well as a putative absence of high en- 
ergy photon absorption on the infrared (IR) back- 
ground ]^^. A review of the various schemes 
can be found in ref. jBU]. There is no firmly es- 
tablished scheme of violating Lorentz invariance 
at present. The most popular scenario invoked in 
explaining the existence of trans-GZK cosmic ray 
events follows the paradigm laid out in Ref. |58j . 
According to that scheme, Lorentz invariance is 
violated by assigning different maximal speeds 
to different particle species (for instance, to elec- 
trons and photons), none of them or some of 
them being equal to what is known to be the 
speed of light in vacuo. 

In this manner, the rate of pair production 
by photons and the rate of photoproduction by 
protons can be suppressed by invoking a rather 
small amount of Lorentz invariance violation. 
Conversely, the absence of trans-GZK flux in 
future experiments could place severe constraints 
on Lorentz invariance violation. It is important 
to point out that the same set of assumptions 
leads to some rather dramatic changes in the 
development of cascade showers, see ref. A 
violation of Lorentz invariance as an explana- 
tion of trans-GZK event cannot be ruled out at 
present. 

4. Conclusion 

A plethora of particle physics models was cre- 
ated to explain the trans-GZK events, which are 
anomalous according to expectations based on as- 
trophysics and on SM. Z-bursts, TD seem to be 
ruled out already to a large extent by the extra- 
galactic diffuse 7-ray background measurements. 

'^We were informed by Michael Fall (private communica- 
tion) that the uncertainties in our knowledge of the IR 
background are considerable. 
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Superheavy dark matter is fatally wounded ac- 
cording to the available observational data on the 
ratio of UHECR photons to protons and on the 
lack of directional asymmetry. Strongly interact- 
ing neutrinos and Lorentz invariance violation are 
still alive. The nature and origin of EHECR is 
still an enigma. If the observational data from 
the next generation of detectors support the ex- 
istence of trans-GZK events then the next great 
leap for particle physics could come from non- 
accelerator experimental physics. It is evident 
that our physics and astrophysics ideas can truly 
be changed by the results coming from the con- 
tinuing and the newly commissioned detectors: 
HiRes, Telescope Array, P. Auger Observatory, 
EUSO, OWL, AMANDA, ICECUBE, NESTOR, 
NEMO, ANTARES, HESS, VERITAS, MAGIC, 
GLUE, FORTE etc. 
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